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摘 要 本文 利 用 CFD 数值 模拟 软件 对 当量 直径 为 200 hm 的 不 同和 孔隙 率 的 又 排 圆 形 微 针 肋 热 沉 的 流动 和 传 热 特 性 进 
行 了 数值 模拟 . 模拟 结果 表明 ， 随 着 雷诺 数 的 增加 ， 针 肋 尾 部 出 现 不 同形 态 的 涡 . 不 同 孔 隙 率 的 通道 中 ,由 于 后 排 针 肋 的 
影响 ， 针 肪 尾 迹 区 的 涡 的 形态 、 大 小 出 现 不 同 。 对 不 同和 孔 际 率 的 微 针 肋 热 沉 ， 热流 密度 的 变化 对 热 沉 的 努 塞 尔 数 Nu 的 
影响 不 同 , 孔隙 率 较 小 时 热流 密度 对 热 沉 的 努 塞 尔 数 Nu 影响 不 大 , 孔隙 率 较 大 时 热流 密度 升 高 , 热 沉 的 努 塞 尔 数 Nw 升 高 . 
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Abstract In this paper, by utilizing CFD, numerical simulations and analyses were used to research 
the flow and heat transfer characteristic in different porosity staggered circular shaped micro pin-fin 
heat sinks whose hydraulic diameter is 200 um. The results indicate that, with the increase of 
Under 
different porosity, the size and shape of vortexes are different due to the influence of pin-fins at 


Reynolds number, vortexes with different structures are found in the pin fins wake zone. 


downstream. The influence on Nusselt number of the heat sinks caused by the change of heat flux is 
different under different porosity. There is barely influence on Nusselt number under smaller porosity 
when heat flux changes, while under larger porosity, the Nusselt number increases with the increase 


of heat flux. 
Key words numerical simulation; micro pin-fin; heat transfer characteristic; low distribution anal- 
ysls 
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随 着 微 电 子 技术 的 快速 发 展 ， 微 流动 系统 在 高 
热流 密度 微型 咒 件 及 设备 的 散热 中 得 到 越 来 越 广泛 
的 应 用 . 作为 微 流动 系统 的 重要 组 成 部 分 , 微 通道 热 
沉 因 具 有 较 好 的 强化 传 热 效 果 ,得 到 了 国内 外 学 者 
的 广泛 关注 和 研究 ， 并 发 展 出 强化 传 热 效 果 更 好 的 
微 针 肋 热 沉 . Kosar 等 0 和 首先 开展 了 微 针 肋 热 沉 
流动 换 热 方 面 的 工作 ， 对 微 针 肋 热 沉 的 摩擦 系数 的 
计算 进行 了 研究 ,提出 了 适用 于 微 针 肋 热 沉 的 摩擦 
系数 关联 式 . 随后 , 用 模拟 和 实验 的 方法 , 对 不 同形 
状 、 不 同 布置 方式 的 微 针 肋 热 沉 的 传 热 与 压 降 进 行 
了 研究 , 指出 通过 合理 地 设计 微 针 肋 热 沉 , 可 以 在 保 
持 优 越 传 热 性 能 的 基础 上 有 效 减 小 流动 阻力 ， 进 而 
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降低 泵 功 消耗 , 周明 正 等 站 和 夏 国 栋 等 @ 分 别 对 流 
体 模 擦 水滴 形 微 针 肋 热 沉 和 获 形 、 长 萎 形 微 针 肋 热 
沉 进行 了 实验 研究 , 研究 表明 ,在 相同 布置 方式 下 ， 
相对 于 贺 形 和 萎 形 针 肪 ， 水 滴 形 和 长 萎 形 针 肋 在 设 
计 上 避免 了 尾 涡 澡 留 区 的 形成 具有 较 好 的 传 热 效 
有 果 。 由 于 在 微 电 子 设备 的 实际 应 用 中 ， 屁 件 的 发 热 
量 会 随 工 作 条 件 的 变化 而 变化 ， 对 热 沉 的 散热 效果 
产生 影响 , 因此 一 些 研究 人 员 也 开始 研究 热流 密度 
变化 对 热 沉 换 热 性 能 的 影响 . Guan 等“- 引 对 不 同 
热流 密度 下 圆 形 针 肋 尾 部 旋涡 的 变化 ， 以 及 热 沉 中 
边界 层 厚度 的 变化 进行 了 研究 ， 发现 热流 密度 升 高 
会 使 针 肋 尾 迹 区 拉 长 ,边界 层 减 薄 ， 壁面 效应 减弱 ， 
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使 热 沉 对 流 换 热能 力 增强 . 

当前 ， 对 微 针 肪 热 沉 的 研究 多 集中 在 传 热 和 压 
降 特 性 方面 ， 对 不 同 排 布 方式 下 多 圆柱 绕 流 的 对 比 
研究 还 较 少 ; 同时， 热流 密度 变化 对 不 同 排 布 方式 
的 微 针 肋 热 沉 的 影响 也 不 明确 。 因 此 ， 本 文 用 数值 
模拟 的 方法 , 在 不 同 热流 密度 下 , 对 当量 直径 为 200 
um, 孔隙 率 为 s=0.743 和 s=0.836 的 两 种 微 针 肋 热 
沉 中 针 肋 扰 流 流 场 和 传 热 特性 进行 了 研究 。 
1 数值 模拟 
1.1 物理 模型 

考虑 到 两 种 微 针 肋 热 沉 模型 为 对 称 结构 ， 因 此 
选取 模型 的 一 半 作 为 模拟 单元 , 模拟 基本 单元 如 图 1 
所 示 , 针 肋 正面 结构 如 图 2 所 示 . 针 肋 排 布 方式 如 表 
1 所 示 . 根据 不 同 的 针 肋 排 布 方式 , 两 个 圆 形 微 针 肋 
热 沉 的 孔 际 率 分 别 为 e。 = 0.743 和 < = 0.836, 将 两 个 
微 针 肪 热 沉 分 别 命名 为 SCP-0.743 和 SCP-0.836( 其 
中 S 代表 又 排 布置 , C 代表 圆 形 , P 代表 针 肋 , 数字 
代表 和 孔 际 率 )。 固体 基板 高 豆 =0.15 mm， 加 热 面 长 
了 =10 mm, 宽 W=1.76 mm, 针 肪 高 H2=200 Um, 前 
后 各 有 3 mm 的 稳 流 区 , 不 进行 加 热 . 


(a) SCP-0.743，s =0.743 
图 1 基本 计算 单元 


Fig. 1 Basic computing elements 


(b) SCP-0.836, se=0.836 


SL 
图 2 微 针 肪 热 沉 正面 示意 图 


Fig. 2 Diagram of the micro-pin fn heat sink 


表 1 针 肋 间距 尺寸 


Table 1 Dimension of pin fin spacing 


模型 ST/mm SL/mm Sp/mm 
SCP-0.743 0.44 0.28 0.36 
SCP-0.836 0.44 0.44 0.49 
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1.2 控制 方程 

为 简化 计算 ,本 文 对 模型 做 出 如 下 假设 : 三 维 、 
稳 态 、 无 内 热源 层 流 流动 和 传 热 过 程 ， 工 质 颖 度 随 
温度 变化 按 分 段 线 性 变化 ; 忽略 辐射 和 热 损失 作用 . 
由 假设 可 知 ， 流 体 的 连续 性 方程 、 动 量 方程 和 能 基 


方程 可 简化 为 ; Dou 


pm (1) 
Ootuiuj) DOP) 0 Ou; 0ui 
Oz; 加 Oz; Dp HE Be | Ba; 局 
O(pruicp#T) 0 oT 
Oz; ~ Oz; 人 外 
(3) 
2 Ou; Ous , Ous 
ES Ozi Oz; Oj 
固体 区 域 能 量 方 程 为 : 
O oT 
re 4 


其 中 , pr 为 流体 密度 ,cp,t 为 流体 比热容 ,Xi 为 流体 
导热 系数 ,Ht 流体 运动 务 度 ，》Xs 为 固体 导热 系数 ， 
1.3 边界 条 件 

假设 流体 与 固体 接触 面 没 有 速度 滑 移 。 边界 条 


件 设 置 为 : 
入 DI: 
T=0, u=un, v=0, w=0, T=Th, (5) 
出 口 
z=L, P=0 (6) 
对 称 边 界 : 
Ou_ o_Ow_ Ou 
Oy Oy Oy Oy 
底面 加 热 面 加 以 恒定 热流 : 
oT 
-hq (8) 
其 他 壁面 均 为 绝热 边界 。 
U= uin, v=0, w=0, a (9) 
On 


2 模拟 结果 与 讨论 
2.1 传 热 特性 
Nu 数 是 表示 对 流 换 热 强度 的 无 量 纲 数 , 可 由 式 


(10) 得 出 : 
_hd 


Nw 一 一 
u 


(10) 
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其 中 , h 为 对 流 换 热 系数 ，W/(n2K)，d 为 当量 直 
径 , m。 

由 于 在 实际 应 用 中 ， 微 型 器 件 和 设备 的 工作 环 
境 、 条 件 和 状态 会 发 生变 化 ， 使 得 散热 器 的 热 负 荷 
发 生变 化 ， 因 此 本 文 研 究 了 不 同 热流 密度 下 微 针 肋 
热 沉 Nu 数 的 变化 情况 . 图 3 所 示 为 SCP-0.743 在 
底面 热流 密度 分 别 为 =4x105 W/m?, 8x105 W/m? 
1.2x106 W/m? 时 ，Nu 数 与 Re 数 变 化 关系 的 曲 
线 图 。 由 图 3 中 可 以 看 出 ， 相 同 热 流 密度 下 ， 随 
着 Re 数 的 增 大 ，Nwu 数 也 逐渐 增 大 ; 相同 Re 数 
下 ， 热 流 密度 的 升 高 ， 对 Nu 数 影响 不 大 . 图 4 所 
示 为 SCP-0.836 在 底面 热流 密度 分 别 为 9=4x105 
W/m?, g=8x105 WA/mm2，1.2x105 W/m? 时 ,Nu 
数 与 Re 数 变 化 关系 的 曲线 图 。 从 图 4 中 可 以 看 
出 , 在 相同 Re 数 下 ， 随 着 热流 密度 升 高 ，Nu 数 有 
较为 明显 的 上 升 。 由 此 可 见 ， 在 相同 Re 数 下 ， 对 


—s—g=4x10” W/m 
—o— g=8x105 W/m 
—A— g=1.2x10 W/m 


40 60 80 100 120 140 


Re 
图 3 SCP-0.743 在 不 同 热 流 密度 下 Nw 数 与 Re 数 关 系 
Fig. 3 Variation of Nu with Re in SCP-0.743 
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—s_4g=4x10” W/m 
—o—g=8x105 W/m 
A g=1.2x10 W/m? 


40 60 80 


100 120 140 
Re 
图 4 SCP-0.836 在 不 同 热 流 密度 下 Nw 数 与 Re 数 关 系 
Fig. 4 Variation of Nu with Re in SCP-0.836 
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不 同 空隙 率 的 微 针 肪 热 沉 ， 热 流 密度 对 热 沉 Nu 数 
的 影响 有 较 大 区 别 . 由 于 流体 流 型 对 热 沉 的 传 热 特 
性 有 很 大 影响 ， 因 此 ， 本 文 将 研究 不 同 孔 际 率 的 微 
针 肪 热 沉 中 流体 流 型 随 热 流 密度 的 变化 ， 来 分 析 流 
型 可 能 对 换 热 产生 的 影响 ， 
2.2 流动 特性 

为 方便 观察 通道 中 流体 的 流 线 图 随 Re 数 和 热 
流 密度 g 升 高 而 变化 的 趋势 ， 本 文 将 通道 沿 流 动 方 
向 ,在 入 口 、 中 段 和 出 口 处 分 别 截取 了 相同 长 度 的 
观察 区 域 . 图 5 为 入 口 雷诺 数 Reis 二 60， 热 流 密 度 
gq 二 4x105W/m? 时 , 热 沉 通道 中 流体 的 流 线 图 . 可 以 
看 出 , 在 Rein=60, 时 , 针 肪 尾部 出 现 一 对 较 小 的 旋 
涡 。 由 于 SCP-0.743 中 后 排 针 肋 距 前 排 针 肋 较 SCP- 
0.836 的 近 , 因此 , 针 肪 尾部 的 尾 涡 区 范围 较 小 , 在 针 
肪 区 出 口 处 , 由 于 SCP-0.743 的 通道 最 大 流速 (1.03 
m/s) 比 SCP-0.836 的 最 大 流速 (0.94 m/s) 大 , SCP- 
0.743 最 后 一 排 针 肋 的 尾 涡 区 明显 长 于 SCP-0.836.， 


3 


N 

(a) SCP 一 0.7 

图 5 Rein = 60, gq = 二 4x105 W/m? 时 , 通道 中 流体 的 流 线 图 
Fig. 5 Flow field in different heat sinks at 

Rein = 60, g = 4 x 105 W/m? 


随 着 Re 数 逐 渐 增 大 ， 针 肋 对 流体 的 扰动 加 剧 ， 
出 现 不 同形 态 的 涡 。 如 图 6 所 示 ， 可 以 观察 到 对 称 
涡 、 单 个 涡 、S 型 涡 等 . 不 同 空 际 率 下 ,相同 观察 区 
域内 的 涡 的 形态 也 不 同 。 

如 图 6(aj)、(b), 当 Reiw=100 时 , 在 SCP-0.743 
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中 可 以 观察 到 少量 单个 涡 、S 型 涡 的 存在 . 而 在 SCP- 
0.836 中 只 有 对 称 涡 存 在 。 随 着 Re 数 的 增 大 ， 如 图 
6(c)、(d)， 当 Rein=160 时 ，SCP-0.743 中 出 现 单 个 
涡 、S 型 涡 的 范围 增加 .。 相同 Re 数 下 , SCP-0.836 也 
开始 出 现 单个 涡 、S 型 涡 , 如 图 6(e)、 (f), 当 Rein=200 
时 , 两 种 热 沉 中 可 以 更 加 明显 观察 到 不 同形 态 的 涡 . 
出 现 这 一 现象 的 原因 是 , 对 于 SCP-0.743, 由 于 孔隙 
率 较 小 , 针 肋 排 布 较 密 , 最 大 速度 较 大 , 扰动 更 加 剧 
烈 , 相对 于 SCP-0.836, 较 小 的 Re 数 下 绕 流 针 肋 产 
生 的 扰动 已 不 能 被 务 性 力 抑制 ， 使 对 称 的 涡 演 变 为 
单个 涡 、S 型 涡 。 随 着 Re 数 的 增 大 , 通道 中 最 大 速 
度 的 增加 使 流体 的 扰动 加 剧 ， 能够 更 多 的 观察 到 针 
肋 尾 迹 区 的 涡 由 对 称 涡 变 为 单个 涡 或 S 型 涡 . 


A 
(e) 
Reis=100, 7 mmL<8.96 mm, (a) SCP—0.743, (b) SCP—0.836; 
Reis=160, 3.64 mm<L<5.6 mm, (c) SCP—0.743, (d) SCP—0.836; 
Rein=200, 3.08 mmL<5.04 mm, (e) SCP—0.743, (f) SCP 一 0.836; 
图 6 不 同 Re 数 下 流 形 的 变化 情况 : 
Fig. 6 Flow field under different Re number: 


由 于 在 实际 应 用 中 散热 右 的 热 负 茶会 发 生变 化 ， 


因此 还 需 研 究 不 同 热 流 密度 下 流体 在 通道 中 流动 时 
流 型 的 变化 。 本 文 用 数值 模拟 的 方法 ,研究 了 相同 
Re 数 时 ， 三 种 不 同 热 流 密度 下 通道 内 流体 流 型 的 
变化 情况 . 图 7 为 SCP-0.743 在 Reis 一 200 时 , 不 
同 热流 密度 下 流 型 的 变化 情况 . 对 比 图 7(a)~(c) 可 
以 看 出 ， 随 着 加 热 面 热流 密度 的 升 高 ， 针 肋 尾 部 涡 
发 展 的 速率 发 生变 化 。 以 出 口 段 为 例 ， 可 以 观察 到 
随 着 热流 密度 的 升 高 ， 尽 管 流体 条 度 减 小 ， 流 体 绕 
流 针 肋 时 黏 清 力 减弱 ， 使 壁面 效应 减弱 ， 可 以 加 强 
换 热 ， 但 由 于 前 后 排 针 肋 间距 较 小 ， 对 流体 形成 挤 
压 , 使 涡 的 出 现 相 对 沾 后 ,流体 扰动 的 剧烈 程度 下 
降 ， 削 弱 了 换 热 能 力 。 两 因素 的 共同 作用 ， 使 热流 
密度 对 Nu 数 的 影响 相对 较 小 . 图 8 为 SCP-0.836 
在 Rein=200 时 ,不同 热流 密度 下 流 型 的 变化 情况 . 
同样 以 出 口 段 为 例 , 可 以 观察 到 随 着 热流 密度 上 升 ， 
流体 温度 上 升 ， 黏 度 下 降 ， 流 体 绕 流 针 肋 时 黏 清 力 
减 小 ， 使 流体 绕 流 针 肪 时 涡 生 成 的 速率 加 快 ， 意味 
着 流体 扰动 的 加 强 , 增强 了 流体 混合 ， 加 强 了 对 流 
换 热 , 因此 Nw 数 随 热 流 密 度 升 高 而 有 相对 明显 的 
升 高 。 


(NY fg WN 人 
A \ / SN _ 
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(a) g=4x105 W/m? (b) g=8x105 W/m?(c) g=1.2x10° W/m? 
图 7 SCP-0.743 热 沉 中 流 型 随 热 流 密 度 变化 情况 : 
Fig. 7 Variation of flow field with heat flux in SCP-0.743 
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图 8 SCP-0.836 热 沉 中 流 型 随 热流 密度 变化 情况 : 
Fig. 8 Variation of flow field with heat fux in SCP-0.836 


3 结 论 


本 文 借助 CFD 模拟 软件 , 用 数值 模拟 的 方法 对 
当量 直径 为 200 wm、 和 孔 际 率 为 e =0.743 和 = =0.836 
Pea etn ee aa 

通过 分 析 模 拟 结 果 得 出 了 以 下 结 

1) 不 同 孔 隙 率 下 ， 0 
能 的 影响 不 同 . 孔隙 率 较 大 时 ， 随 着 热流 密度 上 升 ， 
热 沉 的 Nu 数 有 较为 明显 的 上 升 ; 孔隙 率 较 小 时 , 热 
流 密度 对 热 沉 Nu 数 几乎 没有 影响 ， 

2) Rein=60 时 ， 针 肋 的 尾部 可 以 形成 一 对 对 称 
的 涡 。 当 Re 数 继续 上 升 ， 流 体 绕 流 产生 的 扰动 已 
不 能 被 务 性 力 抑制 ， 在 针 肋 尾部 出 现 对 称 锅 、 单 个 


涡 、S 型 涡 等 不 同形 态 的 涡 。 随 着 Re 数 继续 增 大 ， 
单个 涡 、S 型 涡 出 现 的 更 多 . 孔 际 率 小 的 微 针 肋 热 沉 
中 更 容易 出 现 不 同形 态 的 涡 。 

3) 不 同 孔 际 率 下 ， 热流 密 度 对 流 型 的 影响 有 区 
别 . 孔隙 率 为 0.743 时 , 热流 密度 升 高 会 使 针 肋 尾部 
的 涡 滞 后 形成 ,流体 扰动 减弱 ， 抵 消 了 部 分 流体 秋 
度 降低 、 壁 面 效 应 降低 带 来 的 强化 传 热 效 果 。 孔 际 
率 为 0.836 时 , 随 着 热流 密度 升 高 , 针 肋 尾部 的 涡 加 
速 形成 ， 加 强 了 通道 内 流体 的 混合 和 扰动 ， 带 来 相 
对 较为 明显 的 强化 传 热 效 果 ,Nu 数 有 所 提升 。 
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